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1.  INTRODUCTION 


In  thi  Report  we  investigate  the  scattering  of  scalar  waves  from  a  distribu¬ 
tion  of  a  finite  number  of  scatterers.  The  technique  we  employ  was  developed  by 
Montroll  and  West*  to  describe  the  discrete  scattering  of  scalar  waves  from 
defects  on  regular  lattices.  The  method  makes  no  assumptions  about  the  sym¬ 
metry  of  the  scatterers  and  therefore  can  be  applied  to  inhomogeneities  of  arbi¬ 
trary  shape.  The  majority  of  the  approaches  used  in  the  past  have  had  to  either 

o 

truncate  the  multiple  scattering  integrals  at  some  appropriate  order  or  res- 

Q 

trict  the  investigation  to  scatterers  having  a  high  degree  of  symmetry.  We 
avoid  both  these  limitations  here  by  characterizing  a  scatterer  by  a  distribution 
of  N  defects  and  solving  the  scattering  problem  for  these  N  defects  exactly.  The 
only  limitation  of  the  technique  is  the  maximum  number  of  defects  one  can  use 
to  specify  the  scatterer.  which  in  turn  is  determined  by  limitations  in  computa¬ 
tion  time. 

The  problem  posed  by  the  scattering  of  a  scalar  wave  from  a  fixed  obstacle 
of  known  shape  and  composition  is  primarily  numerical  and  involves  the  solution 
of  a  linear  integral  equation.  The  propagation  of  an  ultrasonic  wave  through  a 
homogeneous,  isotropic  elastic  material  is  described  by  a  scalar  wave  equation. 
The  scattering  of  such  a  wave  from  an  imperfection  in  the  solid  generates  both 
longitudinal  and  transverse  waves,  which  satisfy  linear  integral  equations  with  a 
propagation  kernel  given  by  a  Greens  tensor.^-  At  present  a  numerically  tract¬ 
able  theory  for  the  scattering  of  an  elastic  wave  from  a  flaw  of  arbitrary  shape, 
i.e.  one  not  possessing  spatial  symmetry,  does  not  exist.  In  nondestructive 
evaluation  (NDE),  where  such  waves  are  used  as  probes  to  determine  the  shape, 
size  and  composition  of  flaws  in  materials,  the  analysis  has  been  restricted  to 
either  long  or  short  wave  approximations.  In  the  former  case  the  Bom 


V, 
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approximation  has  been  used  extensively  and  in  the  latter  a  generalization  of 

a 

the  diffractive  geometric  optics  approximation  originally  due  to  Keller  has  been 
used.  An  approach  which  can  span  the  gap  in  wavelengths  between  these  two 
extremes  is  highly  desirable. 

Herein  we  limit  our  discussion  to  scalar  waves  scattering  from  inhomo¬ 
geneities  and  examine  properties  of  the  scattered  wavefield  in  the  feu*  field 
region.  Among  the  scattering  configurations  we  examine  are  line  scatterers  of 
finite  length,  i.e.,  hairline  cracks,  including  the  effects  of  bends  in  the  crack  and 
cracks  of  finite  lateral  extent  also  rectangular  patches  of  scatterers.  The  exact 
scattered  wavefield  from  such  theoretically  repelling,  but  physically  interesting 
objects  has  not  been  calculated  previously. 

In  Section  2  we  review  the  mathematical  formalism  for  the  scattering  of 
scalar  waves  from  defects  on  a  simple  cubic  lattice.  The  lattice  is  chosen  for 
computational  convenience  and  is  not  to  be  confused  with  the  crystal  lattice  of 
the  material.  A  real  metal  such  as  iron  consists  of  an  agglomeration  of  crystals 
having  a  characteristic  dimension  of  0.025  cm  with  arbitrary  shapes  and  orienta¬ 
tions.  Each  of  these  crystals  is  called  a  grain  and  there  exists  a  distribution  of 
grain  sizes  in  such  polycrystalline  materials.  The  distribution  of  grain  sizes 
determines  the  fracture  micro-mechanics  of  the  elastic  material  and  can  dom¬ 
inate  such  effects  as  crack  formation  and  propagation.  Here  we  are  not  con¬ 
cerned  with  the  generation  or  evolution  of  cracks,  but  rather  with  their  detec¬ 
tion.  Thus  we  feel  that  the  use  of  our  computational  lattice  is  justified  for  this 
preliminary  investigation  of  S-wave  scattering  in  NDE.  We  show  in  Section  2  that 
in  the  solution  for  the  scattering  of  a  wave  from  N  defects  no  assumptions  are 
made  about  the  relative  positioning  of  the  defects.  An  extended  scattering 
object  is  therefore  represented  by  a  distribution  of  N  such  defects.  The 
Montroll-West  model  solution  is  not  an  exact  representation  of  the  scattering 


from  a  flaw,  because  the  flaw  is  represented  by  only  N  defects.  However,  the 
ratio  of  determinants  is  the  exact  scattering  solution  to  the  N  point  represen¬ 
tation  of  the  flaw.1  The  exact  nature  of  the  model  enables  one  to  exploit  the 
method  to  determine  the  effects  of  interference  and  resonance  in  the  specular 
reflection  from  cracks. 

In  Section  3  we  discuss  the  physical  interpretation  of  the  form  of  the  solu¬ 
tion  to  the  N  defect  scattering  problem.  One  can  relate  the  terms  in  the  expan¬ 
sion  of  the  determinants  to  scattering  diagrams  and  associate  elements  with 
particular  multiple  scattering  effects.  In  this  way  various  effects  can  be  sys¬ 
tematically  suppressed  by  using  approximations  to  the  expanded  form  of  the 
determinants.  It  is  suggested  that  this  technique  can  be  used  to  check  other 
approximation  methods  in  situations  where  exact  analytic  or  numerical  calcula¬ 
tions  cannot  be  made  using  more  standard  methods. 

In  Section  4  we  present  the  results  of  our  calculations. 
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2.  DISCRETE  SCATTERING  MODEL 
A.  Continuum  Equations 

The  propagation  of  a  scalar  wave  through  a  homogeneous,  isotropic  elastic 
material  is  described  by  the  wave  equation 

V*Kr.t)-±-  -g-  r.t)  =  0  (2.1) 

where  i>{  r.t  )  is  the  wavefield  at  position  r  and  time  t  ;  nd  c  is  the  wave  speed  in 
the  material.  For  a  wavefield  of  frequency  u  ,  i.e.  r.t)  -  $  (  r)  exp  [  - iot  ], 
we  cam  rewrite  (2.1)  as  the  Helmholtz  equation 

~  V"  $  (  r)  +  n2*8*  (  r  )  =  0  (2.2) 

where  the  factor  of  5-has  been  introduced  for  later  convenience  and 

n2  =  o2  /  2czk§  The  Helmholtz  equation  is  an  elliptic  partial  differential  equa¬ 
tion.  so  that  in  order  to  find  the  field  $  at  a  given  point  r  one  must  solve  the 
equation  for  the  field  at  all  points  in  space.  This  property  of  elliptic  partial 
differential  equations  severely  limits  the  analytical  and  numerical  tractability  of 
(2.2). 

The  Helmholtz  equation  arises  in  many  contexts,  e  g.  in  the  propagation  of 
electromagnetic  waves  in  media  where  polarization  effects  are  not  important; 
the  Schrddinger  description  of  a  matter  wave  in  an  energy  eigenstate,  etc.  Each 
of  these  applications  has  precipitated  its  own  sequence  of  approximations  for 
solving  (2.2),  including  the  conversion  of  (2.2)  into  a  parabolic  differential  equa¬ 
tion.  The  physically  interesting  situation  is  when  the  medium  is  non- 
homogeneous,  e.g.  flaws  in  an  elastic  material.  In  this  case  we  write 

n8  =  1  +  K(r)/  kl  (2.3) 

where  V  (  r)  characterizes  the  inhomogeneity.  This  function  can  also  be  the 

potential  in  quantum  mechanics,  the  variable  index  of  refraction  in  a  medium 

with  changing  optical  properties,  the  spatially  dependent  dielectric  coefficient 
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for  a  plasma  with  density  variations,  etc.  The  solution  of  (2.2)  with  (2.3)  inserted 
is  given  by  the  integral  equation 

*  (  r)  =  $o  (  r)  -  f  g{  r,  ^;n0)  V  ( t )  $  (  r1)  dV  (2.4) 

where  g  (  r,  r';n0)  is  the  Greens  function  that  satisfies  the  equation 

^g{  r.  r'mo)  +  k§  g  (  r,  r'\n0)  -  6  (  r  -  r'),  (2.5) 

and  $o(  r)  is  the  wavefield  when  V(  r )  =  0. 

B.  Discrete  Equations 

The  discrete  analog  of  the  Helmholtz  equation  is  here  obtained  by  interpret¬ 
ing  the  V*  in  (2. 1)  as  the  second  difference  "Laplacian" 

*  s  ~T  t  (Ej-2  +  Ef1)  (2.6) 

o  >=t 

where  the  Efvs  are  the  discrete  translation  operators  on  a  simple  cubic  lattice, 
i.e. 

f  (lv  li  ■  In)  »  •  Li-1-  In) 

=  Efxnix---li  +  1  -  In)-  (2.7) 

The  position  vector  r  =  1 6  has  the  lattice  components  (I1.i2.i3)  with  the  ele¬ 
ments  of  i  being  integers  and  b  is  the  lattice  spacing.  The  free  lattice  can 
support  the  propagation  of  plane  waves  with  frequency  u  and  wave  vector 
k  ss  (&i,Jfc2,fca)  so  we  substitute  the  plane  wave 

$o(l  *)  -  expji  k  '  I6(  (2.8) 

into  (2.2)  and  replace  nzk§  by  the  eigenvalues  X*.  i.e. 

X*  =  pr[  3  -  cos(bfci)  -  cos(bfc2)  -  cos(bka)].  (2.9) 

We  impose  periodic  boundary  conditions  on  the  lattice  such  that 
$  (*  +  Nb.y.z)  =  $(x,y  +  Nb,z)  =  ${z,y,z  +  Nb)  =  #(x,y,z)  .  This  imposi¬ 
tion  restricts  the  components  of  k  to  integral  multiples  of  2-n  /  Nb,  i.e., 

2n 

k  =  -r^-n  where  n  =  (ni,n2,na)  and  \rii.nz,n^\  =  0,1,  •  •  •  N.  Introducing  the 
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lattice  Greens  function 

G  ( 1  6  ;n0)  =  Am  exp \i  ^-n  •  1{  (2.10) 


into  the  discrete  analog  of  (2.5)  and  equating  coefficients  of  exp  \i  2tr  n  .  1  /  N  j, 
yields  A  a  =  (n,2  -  X*)-1.  Thus  the  Greens  function  on  the  simple  cubic  lattice  is 

given  by 


G  (  U»  .n0)  =  £ 

n  n,.n8.n3  =  l 


jk02  62 ' 3) 


exp  H 2tt  nl/N\ _ 

‘  3]  +  £  cos  (2  nrij/  N) 

1  1 


(2.11) 


when  there  are  N  lattice  sites  in  each  of  the  three  dimensions.  The  Greens  func¬ 
tion  (2.11)  describes  the  free  propagation  of  a  plane  wave  on  a  homogeneous  iso¬ 
tropic  simple  cubic  lattice. 


C.  "Potential"  Scattering 

We  consider  a  set  of  time  independent  defects  on  the  set  of  lattice  points 
a  s  1*  2  ....  n.  In  the  presence  of  these  defects  the  discrete  Helmholtz 
equation  takes  on  the  inhomogeneous  form 

|-Az  +  fc£  *  (16  )  =  -XK(1)$(16)  (2.12) 

where 

,  x  l  =  U.  i  , 

^  ,  0  otherwise  ^ 

and  X  (  s  is  the  strength  of  the  defect,  here  produced  by  a  variation  in  the 
Po 

density  at  the  set  of  sites  { 1  «  }.  We  shall  refer  to  V  ( 1 )  as  the  scattering  poten¬ 
tial  for  the  wave  and  X  as  its  strength.  The  plane  wave  (2.8)  is  a  solution  to  (2.12) 
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in  the  limit  X  -*  0.  The  full  solution  of  the  inhomogeneous  equation  (2.12)  is  then 
a  linear  superposition  of  the  incident  plane  wave  and  the  scattered  wave  pro¬ 
duced  by  the  potential,  i.e. 

$(16)  =  exp{ikj-16J  -  x£  G(l-l„;n0)  K(l«)*(l«6)  (2.14) 

«  =  l 

which  is  analogous  to  the  continuum  solution  (2.4).  We  have  indicated  the  value 
of  the  incident  wave  vector  by  k  *  in  (2. 14). 

To  solve  equation  (2.14)  for  the  total  wave  field  $  (1  6)  we  introduce  the 
notation 

Gaft  =  C(la-lp  -,  n0)  ;  Coa  s  G(  1  - 1  a  ;n0)  ;  $a  =  $( la6 ) 

aliu 

V«sV(Ia).  (2.15) 

Then  if  we  successively  let  1  be  equal  to  the  location  of  each  of  the  n  potential 

sites  { la  j  in  (2. 14)  we  obtain  the  n  inhomogeneous  linear  equations 
(1  +  XGU  Vi)  $i  +  X  Giz  V2$z . +  X  Gin  Vn  $2  =  exp  jiki  l^b  j 

X  Gzi  +  (l  +^X  Gzz  72)$2 . +  X  Gzn  In  =  exp  Jikj  •  ^6  j 

$i  +  $a- •  •  +  (l+*Qm  Jn)$n  =  exp{iki  I„6  | 

The  solution  to  this  set  of  equations  when  n=2  is 


1 

exp(i  kt  ■  1 

i  6)  X  G\z  72  ] 

(2.17a) 

a2 

exp(ik(  ■  lj 

i  b  )  1  +  } 

<  Ggz  V2  j 

1 

1  +  X  Gji  7, 

exp(  i  k* 

1 1  6  ) 

=  *z 

X  Gz\  V\ 

exp( t  k i  ■ 

lz  b) 

(2.17b) 

where 
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1  +  X  G,,  Vj  X  G j2  Vjj 
X  Ggj  V\  1  +  X  C22  Vg 


(2.17c) 


When  (2.17a)  and  (2.17b)  are  introduced  into  the  solution  (2.14),  the  full  expres¬ 
sion  for  $(  lb)  becomes  a  sum  of  terms  which  is  easily  arranged  as  the  ratio  of 
two  determinants.  The  total  wavefield  at  the  observation  site  1  is  then  given  by 


exp(  i  k  i  •  lb)  XG0|K,  X  Goa  ^ 
$(lb)=  —  exp(  i  k  t  -  1 1  b )  1  +  \GnVr  X  G,z  V2 

exp(  i  k  i  l2b)  X  C2|  1  +  XG22  Vz 


(2.18) 


The  above  argument  may  be  repeated  in  every  detail  to  obtain  the  total 
wavefield  after  scattering  from  n  point  potentials: 

exp(iki  lb)  X  Go,  K,  •  •  •  X  Gon 

.  exp(iki  -  1,  b)  l+XGi,V,  XGlnVn 
$(lb)=^-  j  j  (2.19) 

exp(iki  1„  b)  XC^jK,  1  +  X  G^n  Vn 

where 

^  b  detl  6 op  +  X  G,,#  Vp  I ;  a,/3  =  1,2,  n  .  (2  20) 

We  define  the  function  G,  (1,2,  n)  by  the  determinant 

exp( tkj  lib)  A  G,z  V2  ■  ■  •  X  G,n  Vn 
£>,(  1,2,  -  n)  *  j  1  +  X  G22  V2  :  (2.21) 

exp(  i  kj  ■  lnb  )  ■  1  + X  Ga  Vn 

and  are  therefore  able  to  express  the  wavefield  at  the  site  lj  as 

*(l»b)  =  /?,  (1.2.  n)/An  .  (2.22) 

At  the  arbitrary  potential  site  1  we  can  write  the  wavefield  as 

(2.23) 

where  the  ordering  of  the  indices  in  the  determinant  Dp  is  just  a  cyclic  permuta¬ 
tion  of  the  lattice  site  indices  in  (2.21).  The  exact  total  wavefield  can  then  be 


written  as 


(2.24) 


$(16)  =  expjikilbj  -  ■£—  £  GQaVaDa 

“n  0=1 

which  is  an  expansion  of  the  determinant  solution  (2.19)  in  terms  of  minors. 


3.  INTERPRETATION  OF  EXACT  SOLUTION 

In  this  section  we  investigate  expressions  for  the  scattered  wavefield  as 
sums  of  contributions  from  single  scattering,  double  scattering,  triple  scattering 
etc.  The  series  we  obtain  are  analogous  to  the  Born  series,  but  the  individual 
propagators  are  renormalized,  i.e.  they  contain  factors  that  correspond  to  the 
resummation  of  interaction  terms  of  a  particular  form.  Let  us  first  examine 
various  choices  for  these  resummations. 

A.  Renormalized  Potentials  and  Propagators 

Consider  the  following  factored  form  for  the  determinant  A,,: 

a  (1  +  X  Gn  Fj)  (1  +  X  Ga2  Vj>)  '  •  An  (3.1) 

where  we  define  the  scaled  determinant 

An  a  detldn,  +  (1  -«<*)  X  Gap  Vp  I  ;  a.p  =  1.  2.  •  •  •  n  (3.2) 

and  Vp  is  to  be  interpreted  as  the  renormalized  potential 

Vp  =  Vp{l  +  \GppVp)-'  (3.3) 

«•» 

In  terms  of  An  and  Vp  the  wavefield  (2. 19)  can  be  expressed  as 

exp(ik<  16)  X<?01Fi  X  Gbn  Vn 

,  exp(  i  ki  •  1,  6)  1  XG|n  Vn 

*(i b )  =  ~  :  :  (3.4) 

An  :  : 

exp(  i  k«  -  1*6)  \GniVi  1 

Each  off-diagonal  element  of  the  determinant  (3.4)  has  the  form 

IGtpVp  =  =  KG'pVp-K'G.p  Vp  Gpp  Vp  +  0(\S)  .  (3.5) 

which  corresponds  to  a  self-interaction  term  in  a  diagrammatic  expansion. 

Such  resufmmattons  are  mandatory  when  considering  scattering  from  potentials 


^."l  1LT  r.  «'.  ^ '  ~  ’.V.V.’i*  %-  \* 
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whose  strength  becomes  singular  at  short  distances,  e.g.  nuclear  scattering, 
hard  spheres  collisions  or  specular  reflections  of  scalar  waves  from  caustics.  It  is 
clear  that  as  X  Vp  -* «  each  terms  on  the  right  hand  side  of  (3.5)  diverges, 
whereas  the  infinite  series  converges  to  GaP  /  Gpp  .  Thus  all  the  terms  in  the 
wavefleld  (3.4)  are  expressed  in  terms  of  the  renormalized  potentials  so  that 

each  term  individually  remains  finite  as  \  ,  i.e.  lim  X  Vp  =  1/  Gpp. 

X  Vp 

In  the  second  renormalized  form  the  determinant  An  is  written  in  the  fac¬ 
tored  form 


A„  =  ( 1  +  X  Gn  Vj)  (1  +  X  Ggj V2)  ■  •  •  An 


where  we  define  the  scaled  determinant  A  n  as 


A„  s  det!da,  +  (1  -Stf)  X  Gap  Vp\  ,a,p  =  1.2.  •  •  ■  n  (3.7) 

and  Gafi  is  to  be  interpreted  as  the  renormalized  propagator 

~G«  =  GaP{l+\G9aVa)*(l  +  \GppVp)-*.  (3.8) 

Again  we  note  that  X  tVp  remains  finite  as  the  set  of  potentials  {  Va  j  become 
infinite  (as  long  as  they  all  become  infinite  in  the  same  way). 

*f  ~ 

In  terms  of  A  n  and  GaP  the  total  wavefleld  can  be  written  as 


exp(ikjlb)  Ux 


$(16)  = 


A  C,«  C4 


XGn,£/, 


where  we  have  introduced  the  functions. 


(3.10) 


Ua  =  A  Goa  Va  /  (1  +AGaa  Va)*. 


(3.11a) 


=  exp  (  i  ki  .  1  a  6  )/  (1  +  X  Gm  Va)*. 


(3.11b) 


In  this  factorization,  an  expansion  of  the  determinant  (3.10)  in  minors  always 
yields  terms  of  the  form  Ua  W„  which  remain  finite  for  singular  potentials.  The 
renormalization  Greens  functions  in  (3. 10)  can  be  interpreted  as  the  propaga¬ 
tion  of  a  plane  wave  on  a  lattice  distorted  by  the  presence  of  the 
n  point-potentials. 


E  Multiple  Scattering  Approximations 

Consider  first  the  case  of  plane  wave  scattering  from  three  point  potentials, 
J  =  I1.I2.I3-  The  solution  (3.4)  can  be  expressed  in  terms  of: 


as 


9ap  =  ^  Ogf  Vf 


§(16)  =  ~ 

Aa 


exp(ik<  .16  ) 
exp(ik* .  lj  6 ) 
exp(ik{.l2b) 
exp(ikt  I36) 


9  01 

9  02 

9  03 

1 

9 12 

9  13 

921 

1 

023 

9  31 

9  32 

1 

Expanding  the  determinant  in  the  numerator  of  (3.13)  yields 


(3.12) 


(3.13) 


§(16)  =  e 


-  1* 


90a  ® 


<W 


a=l 


Az(a)  +  £  <7,*  + 

0t»a  0*7 


(3.14) 


The  quantity  A2  {!)  is  the  2x2  determinant  A2  which  involves  only  the  g’s  with 
subscripts  2  and  3  but  not  1.  Each  term  in  the  expansion  can  be  identified  with  a 
scattering  diagram,  [cf.  Figure  l]  .  Figure  la  represents  the  single  scattering  of 
our  original  plane  wave  by  the  potential  at  site  1  v  The  full  form  of  the  contribu¬ 
tion  of  that  term  to  §  ( 16 )  is 


A2(l)  _<k4  1,6 
— - e 

A3 


A  G01 


V\  e 


>1* 


A2(T) 

A3 


(315) 


The  factor  exp(i  k,  li  6)  in  (3.15)  represents  the  plane  wave  incident  upon  the 
scattering  point  1,;  \VX  represents  the  effective  (renormalized)  scattering 


potential  and  G0i  corresponds  to  the  scattered  wave  which  propagates  from  l! 
to  the  observer  at  1.  The  ratio  of  the  A's  is  the  contribution  of  the  other  scatter¬ 
ing  points  to  the  single  scattered  amplitude  even  though  they  are  not  directly 
involved  in  the  scattering.  They  give  rise  to  the  distortion  of  the  underling  lat¬ 
tice  very  much  like  a  mean  field  in  which  the  single  scattering  occurs.  As  X  -* 

•M 

0.  Vx  -»  Vx  [cf.  (3.3)3  and  the  ratio  of  the  A's  approaches  unity  so  that  (3.15) 
becomes  a  contribution  to  the  first  Born  approximation  of  the  solution  in  this 
limit. 

In  Figure  lb  the  diagram  corresponding  to  the  double  scattering  term 


As-1 003  0ai  expfiki-libj  (3.16) 

is  depicted.  Interpreting  the  factors  in  (3.16)  from  right  to  left,  we  have,  the 
plane  wave  incident  on  lattice  site  lj  ,  the  scattering  from  1  j  to  13  followed  by 
the  scattering  from  la  to  the  observer  at  I.  The  influence  of  those  points  not 

directly  involved  in  the  scattering  is  A  i  ( 1, 3)/  A  3.  but  since  A ,  =  1 .  only  the  A  3 
factor  appears. 

These  results  are  immediately  generalized  to  the  case  of  m—th  order 
scattering  in  a  system  of  n  scattering  centers.  The  contribution  to  $  ( lb  )  of 
the  scattering  sequence  li  -» lam_,  -*  •  •  laj  -*  observer  is 


0Oa,0aiaz  0  am_2'Om_i0  a 


<k4 


m-t.l 


1,8  A, 


-rn  •••  .am-1.1)/ An  (3.1?) 


where  we  let  A  n_m  (ai,a8,  •  •  •  ,nm_,,l)  represent  the  determinant  An_m  which 
involves  the  (n-m)  points  not  including  the  points  from  which  the  wave  has  been 
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scattered.  Hence  the  field  observed  at  site  1  can  be  expressed  in  the  series 
form 

$(16)  =  expji  ki  - 16  j  +  £  go*  e'k<  A„_,  (a)/ An 

asl 

+  £  9o*9afi  e<ki  An_2(5.0)  /  An  +  •  (3.10) 

ajl=\ 

This  expansion  is  the  generalization  of  the  multiple  Born  expansion,  but  with  all 
terms  renormalized. 


4.  SCATTERING  CROSS  SECTION 


As  mentioned  in  the  Introduction,  the  use  of  ultrasonic  waves  as  a  probe 
into  the  structure  of  materials  has  been  very  fruitful.  In  particular,  measure¬ 
ments  of  the  total  power  transmitted  from  source  to  receiver  through  an  elastic 
material  provides  information  on  the  distribution  of  scattering  centers  within 
that  material.  The  attenuation  of  the  ultrasonic  wave  in  a  distance  z  along  its 
direction  of  propagation  is  exp  [  -  cr(cJ)z  ]  where  the  attenuation  coefficient  o(o) 
is  dependent  on  the  frequency  of  the  incident  wave  and  is  proportional  to  the 
average  cross  section  of  the  scatterers  in  the  material.  This  relation  has  been 
used  to  determine  the  distribution  of  grain  sizes  in  polycrystalline  materials.  If 
the  attenuation  of  the  wave  is  due  primarily  to  the  elastic  scattering  of  energy 
out  of  the  direction  of  propagation,  then  the  elastic  scattering  cross  section  is  a 
quantitative  measure  of  how  much  the  incident  wave  is  attenuated  along  the 
straight  line  path  joining  the  source  and  receiver.  Even  in  the  case  of  a  localized 
scatterer  such  as  a  crack,  the  scattering  cross  section  determines  the  intensity 
of  the  signal  received  at  different  observation  points. 

In  the  continuum  the  total  wave  field  can  be  written  as 

*(r)  =  expjiki  rj  +  e*Rj?*Z±/  (kt.  k,)  (4.1) 

where  the  incident  wave  vector  is  kt  =  k  Jt0  .  the  final  wave  vector  is  ky  =  rfc0 

A  A 

in  terms  of  the  unit  vectors  k  =  k/  fc ,  r  =  r /  r.  Here  r  is  the  asymptotic  loca¬ 
tion  of  the  observer;  /  ( k4,k/  )  is  the  amplitude  of  the  scattered  wave  and 

cos  -*(  k4  k j  )  =  9  is  the  scattering  angle.  The  differential  scattering  cross 
section  is  given  by 

=  1/ (kt.kj)!2  (4.2) 

and  the  total  scattering  cross  section  is  obtained  by  integrating  (4.2)  over  the 


solid  angle  dfl  =  sin0d0d^  i.c 


atot  =  2-nf  l/(ki,k/)lz  sin  0  d  0 
0 


The  total  wavefleld  on  the  lattice  is  given  by  (2.24)  as 

#(16)  =  expliki-lbi-^S  Goa  Va-Da  (4-4 

The  scattered  intensity  is  observed  at  1  which  is  generally  a  considerable  dis¬ 
tance  from  a  typical  scattering  site  la  .  Hence,  the  large  1  approximation  can  be 
used  for  Go*  .  This  is  just  the  Greens  function  for  the  original  continuum 
medium.  The  other  G's  in  the  determinant  Da  are  given  by  the  appropriate  lat¬ 
tice  Greens  functions  [cf.  (2.11)]  determining  the  wave  propagation  among 
potential  sites  { lfl  J.  We  consider  the  limit  of  the  lattice  Greens  function  (2.11) 

as  N  -*  “  and  define  the  continuous  variable  lim  -rr-  =  <pj  so  that  (2. 1 1) 

TV 

becomes 


explil 


(fc|  6z-3)  +  ^  cos <pj 
i=i 


Now  writing  the  continuum  position  r  =  16  we  can  rewrite  the  lattice  Greens 


function  as 


G(1&;toq)  =  63G0(rfc0)  • 


The  asymptotic  form  of  the  continuum  Greens  function  is 


so  that  from  (4.4)  we  obtain 
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Recall  that  k/  =  k0  r  so  that  the  amplitude  of  the  scattered  wavefield  is 

f(kiif)  =  e"kor  r“  VaDa  (4.9) 

a=l 

since  r  » !  rfll .  The  dependence  of  the  scattering  amplitude  on  k<  is  con¬ 
tained  in  Da- 

In  a  subsequent  section  we  calculate  /(ki.k/)  for  various  numbers  and 
configurations  of  point  scatterers. 
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5.  CALCULATION  OF  LATTICE  GREENS  FUNCTION 

It  is  clear  from  all  our  discussion  that  the  scattering  of  waves  off  defects 
specified  by  Va  on  a  lattice  is  completely  determined  by  the  lattice  Greens  func¬ 
tion.  On  a  simple  cubic  lattice  we  rewrite  (4.5)  as 


*(0;I.m,n)  =  4-  /  dx  /  f  dz 

rr  n  Q-cosx -cosy -cose 


where  the  triad  (i.m.n)  specify  the  integer  lattice  coordinates,  and  G  =  -b  <S. 
Because  the  integrand  of  (5.1)  is  both  singular  and  oscillatory  the  numerical 
integration  is  non-trivial  and  requires  a  certain  amount  of  discussion.  There  is  a 
substantial  literature  on  Greens  functions  of  the  form  (5.1)  since  they  arise  in 
the  model  studies  of  many  physical  phenomena  e.g.  in  the  study  of  the  distribu¬ 
tion  of  magnetization  around  an  impurity,  i.e.  localized  spin  wave  modes  and 
other  problems  in  condensed  matter  physics. 


The  emphasis  in  the  literature  published  on  the  numerical  solutions  of  (5.1) 
♦ 

has  concentrated  on  obtaining  accurate  values  of  ®(0; 0,0,0),  i.e.  the  Greens 
function  at  the  origin.  For  the  purposes  of  scattering  theory  we  need  to  do 
much  better  than  this  and  obtain  accurate  values  of  for  arbitrary  (£,m,n)  . 
Accuracy  is  quite  important  in  order  to  faithfully  calculate  the  effect  of  interfer¬ 
ence  among  the  scattered  wave  components  when  many  defects  are  present. 
Recall  that  the  exact  results  can  be  expressed  as  an  infinite  perturbation  series 
so  that  a  modest  error  in  phase  in  any  iterated  component  could  yield  substan¬ 
tial  spurious  effects. 

An  iterative  technique  has  recently  been  developed  by  Morita  in  which  a 
knowledge  of  the  Greens  function  at  the  three  sites  (0,  0,  0),  (2,  0,  0)  and  (3,  0,  0) 
can  be  used  to  generate  values  of  ®(Q;l  ,m,n)  at  any  of  the  other  lattice  coordi¬ 
nates.  We  do  not  present  Morita's  complete  arguments  here,  but  rather  quote 
the  equations  necessary  to  calculate  the  Greens  function.  In  particular  we  point 


out  an  error  in  one  of  the  puousned  expressions. 

The  recurrence  relations  that  are  used  to  connect  the  value  of  the  Greens 
function  at  the  sites  (0.  0,  0),  (2.  0,  0)  and  (3,  0.  0)  to  its  value  at  other  lattice 
sites  is  most  easily  written  in  a  rotated  system  of  coordinates.  We  rotate  the 
angles  x  and  y  in  (4.1)  to  x’  and  y’  as  follows: 

x'  =  (x+y)/2  ;  y"  =  (x -y)/2 

£*  =  l  +  to  ;  ml  =  l  -m  (5.2) 

so  that  the  Greens  function  on  the  rotated  system  are 

^‘(Oii'.m'.n)  =  ^(Q;  ( V  +  nx')/2.(f -m')/2.n)  (5.3a) 

*(Q;  l  +m,l-mlTi)  =^(Q;£.m.n)  (5.3b) 

as  can  be  verified  by  direct  substitution  into  the  integral  (5. 1).  Note  that  V  and 

to'  are  either  both  even,  or  both  odd  integers. 

Moritas  scheme  is  based  on  the  difference  equation  for  the  lattice  Greens 
function  restricted  to  the  n=0  plane.  The  specific  equation  he  obtains  is 

(fix'  +  1  +2,  to'  +2. 0)  +  23?1  (f  .fix’  +  2,0)  +  <GX (r  -2, fix’  +2. 0)] 

— 2(2nx'  +  1)  fl  \S\V  +  l,m'+l,0)  +  «‘(i’  -  l.nx’  +  1,0)] 

+  2 m*  \<g\V  +  2, fix'  ,0)  +  2CJ2  &(r  .fix'  ,0)  +  <s\v  -  2, m.'  .0)] 

-2(2nx'-l)  n  \g\V  +  i.m'— 1,0)  +  <Sx{i  -1, fix’ -1.0)] 

+  (m'-l)[«rM(r  +  2, nx' -2,0)  +  2<S,(r, nx' -2,0)  +  »!(r -2,m' -2.0)]  =  0  .(5.4) 
The  manipulation  of  (5.4)  to  create  a  recursion  formula  for  (GX  can  be  found  in 

Moritas'  paper.  Figure  (2)  depicts  a  flow  chart  of  the  logic  used  in  the  computer 
program  to  obtain  (8X  .  The  quantities  A  through  F  used  in  this  procedure  are: 


A  =  p~l  |  (4p-l)  n  [«l(2p,2,0)  +  «?*(2p.0,0)  j 

' (2p -1)  [«/l(?p-1.3,0)  +  (20*  +  1)  Sl  (2p— 1,1,0)  ] 

+  (4p-3)  n  [« *(2p -2,2,0)  +  <8x(2p -2,0,0)  ] 

-(p-1)  [s‘(2p  -3,3,0)  +  3 <SX  (2p— 3,1,0)] 

B  =  Q[3(8l(2p.2,0)  +  «51(2p,0,0)] 

-2  [<6x{Zp -1,3.0)  +  0 *<8X  (2p -1.1,0)] 

+  n[3«,(2p-2.2.0)  +  '8\2p-2, 0. 0)] 
-[®1(2p-3,3,0)-^1(2p-3.1.0)]  . 

C  =  (m‘  +  l)[29l(<'tm‘  +  2,0)+tf1(r-2,m,  +  2.0)] 

-2n(2m'  +  1)  [«*(£’  +  l.m*  +  1.0)  +  z?1  (£’— l.m'  r  1,  0)] 

+  2m'  [«>(£’  +  2,m’, 0)  +  2fi2^,(f\m',0)  +  «Sl(£'-2.m\0)] 

-2Q  (2m1  -1)  [*»  (i‘  +  l.m’-l.O)  +  (I'-l.m'-l.O)] 

+  (m'-l)  [«l(£'+2,m,-2,0)  +  2«1(£',m’-2,0)] 


0  =  (2p  +  l)-‘ 


2(4p  +  l)n®1(2p-t-l,  1,0) 


-  4p[®1(2p,2,0)  +  n Z(8X  (25o,0,0)  ] 

+  2(4p-l)fi  <8X  (8p -1,1,0) 

-  (2 P -1)  [f8l{2p -2, 2, 0)  +  <3 1  (2 p  -2,0, 0)] 


£  =  2fl[f  #‘(2p  +  l,3,0)  +  3^‘(2p +  1.1.0)] 


-  6  'S'  (2p ,  4,  0)  -  Bf)2  (2p ,  2, 0)  -  2  (2p ,  0. 0) 


+  20[5  (Sl  (2p— 1, 3. 0)  +  32?*(2p-1.1.0)] 


-3»1(2p-2.4,0)-4«I(2p— 2.2.0)  -  ^1(2p-2,  0. 0)  (5.J 

F  =  20[(4+p)tf1(2p  + 1.3,0)  +  (2-p)  <gl  (2p  + 1,1,0)] 

-3^(2?. 4.0)  -4(fl*(p+l)+l)  <gl{2p,2,0)  +  (4p  D*-l)  »l(8p,0,0) 

+  2fl[(l+p)£1(2p -1.3,0)  +  (1-p)  <3X  (2p— 1. 1.0)]  (5. 

Equations  (5.5a)  -  (5.5c)  inclusive,  correspond  to  Eqs.  (5.7).  (5.4)  (5.9b)  and 

(5.10b)  of  Morita  respectively.  However,  Moritas'  Eq.  (5.12b)  is  incorrect  and 

(5.5e)  takes  its  place. 

Finally,  to  calculate  values  of  <&1(l',m\n)  for  nonzero  n  (and  hence 
<8(l.m,n)  also  ).  we  use  a  modification  of  equations  (5.3)and  (5.4)  of  Morita 
namely 

M)  =  n«‘(r.  m’.O)-  +l.m*  +  1.0) 

+  <8\l’  +  l.m'-  1.0)  +  <g'  (f’-i.m’  +  1.0) 

-jg-tfM*'- l.W-1.0)  (5 

and.  for  n  S:  2, 

<gx(l\m',n )  =  20  «1(l,tm\»-l)-*I(l\m\n-2) 

-«*(r  +  l.m’+l.n-l)-  ^‘(r-l.m’-l.n-l) 

“ (i’+l.m'— l.n— 1)  -®1  (f— l.m'-f l.n-l)  (5 

Before  the  above  recursion  scheme  can  be  applied,  however,  we  require 

8(0.  0.  0).  8(2,  0.  0)  and  <8(3,  0.  0)  (  <3l  (0.  0.  0),  81  (2,  2.  0)  and  81  (3.  3.  C 

as  input. 

The  simple  cubic  lattice  Greens  function  <8( 0  ;  i,  0,  0)  has  a  convenient 
representation  in  terms  of  complete  elliptic  integrals  of  the  first  kind.  The 
Greens  function  is  pure  real  for  3  ^  0  .  and  is  given  by 

i  n 

«Z(0&3;i,0,0)  =  -~r  f  dQ  cos  (IQ)kK(k)  (5 
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where  «s 2(f)  -  cos9)_1  is  the  parameter  of  the  complete  elliptic  integral 


rt  Ok 

(s-9> 

o  vl  -  arsur  0 

For  1  ^  0  ^  3,  one  has  for  the  real  and  imaginary  parts  of  the  Greens  function 

cos-'tn-a) 

^*(fi;J,0,0)  =  ~  f  dQcos(lQ)  K  (Kt) 

1*  0 


+  -J-  f  d9cos(£0)A'(/c)  (5.10) 

"  cos-^n-g) 


eos"’(n-2) 

(fiii.o.o)  =  -3-  /  d9  cos(i0)  K (k'i)  (5.11) 

Tr  0 

where 

Ki  =  k~1  and  K\  =  Vl  - Ki  .  (5.12) 

In  the  range  0  ^  fl  ^  i  ,  we  have 

cos_1n  it 

«*(n;I,0,0)  =  J  d@cos{lQ)K(K.  )  r  /  d0cos(f  0)  K(Kt)  (5.13) 

^  0  ^  cos-’n 


1  n 

(n;f,0,0)  =  —f  dQ  cos(i0)  K  (k\)  (5.14) 

rr  0 

where  /Cj  and  are  given  in  (5.12).  Formulas  (5.B)  -  (5.14)  were  numerically 
evaluated  for  l  =  0,  2  and  3  and  used  as  input  for  the  recursion  scheme  previ¬ 
ously  described.  Figure  3  shows  plots  of  the  reed  and  imaginary  parts  of 
«(0  ;0.  0.  0),  fG{  0;  2,  0.  0)  and  (G{Q  ;  3,  0.  0). 


6.  RESULTS  OF  SCATTERING  CALCULATIONS 


In  this  section  we  present  some  results  of  scattering  plane  waves  off  various 
planar  configurations  of  point  scatterers.  We  have  chosen  the  incident  radiation 

to  have  wavelength  equal  to  the  lattice  spacing  b.  and  to  propagate  in  the  z 
direction.  All  scatterers  lie  in  the  perpendicular  x-y  plane.  Polar  angles  0,  tp 
define  the  directions  along  which  we  view  the  scattering  amplitude  f(0,  <p)  (see 
figure  4). 

A  Scattering  From  a  Line. 

The  first  sequence  of  numerical  results  we  show  are  polar  perspective  plots 
of  the  magnitude  I  f  (0,jfl) !  for  one,  two,  three,  four  and  five  scatterers  arranged 
to  form  a  line.  If  just  a  single  scattering  center  is  present  (located  at  the  origin 
(0,  0.  0).  say),  there  are  no  preferred  directions  for  the  scattering,  therefore  If! 
is  uniform  over  all  angles  0,p  .  The  result  of  adding  a  second  scattered  at  (1,  0, 

0)  is  shown  in  Figure  5a.  (The  number  and  position  «f  ti*e  scattfws  is  indicated 
In  the  upper  left  hand  comer  of  the  figure.)  We  see  that  the  scattering  ampli¬ 
tude  in  maximal  in  all  directions  perpendicular  to  the  preferred  axis,  i.e.  along 
=  rr/2  and  3n/  2  for  all  0  values.  I  f  I  is  also  maximal  parallel  to  the  x-axis, 
at  tp  =  0,  7r  and  0  =  n  /  2  .  The  amplitude  decreases  away  from  these  direc¬ 
tions  due  to  destructive  interference  among  ray  paths.  The  addition  of  a  third 
scatterer,  at  (2,  0,  0),  gives  rise  to  Figure  5b.  The  directions  of  maximall  f  I 
remains  the  same  (serving  to  identify  the  line  of  scatterers),  but  I  f  I  decreases 
more  rapidly  away  from  these  regions  than  for  the  case  of  two  scatterers.  More 
significantly,  a  subsidiary  ridge  has  appeared  as  a  new  topographic  feature, 
marking  the  influence  of  multiple  scattering  events.  Figures  5c  and  5d  show  the 
effect  of  adding  a  fourth  and  a  fifth  defect.  We  see  a  further  sharpening  of  the 
main  ridges,  and  the  introduction  of  additional  subsidiary  ridges. 
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B.  Scattering  from  Perpendicular  lines 

Figures  6a  -  6d  show  a  sequence  of  plots  of  I  f  (  Q,tp  )  i  as  we  add  to  the  line, 
but  now  in  a  perpendicular  direction.  This  sequence  is  most  easily  understood  if 
we  first  consider  the  scattering  off  a  line  of  defects  arranged  along  the  y-axis 
(instead  of  the  arrangement  along  the  x-axis  which  led  to  Figure  5d).  Symmetry 
implies  that  such  an  arrangement  should  lead  to  a  scattering  amplitude  f  (  0,  ip  ) 
related  to  the  previous  amplitude  by  the  rotation  0.  <p  -*  Q.<p  +  tt/  2.  All  features 
of  the  perspective  plot  of  I  f !  would  simply  be  shifted  by  n/  2  along  the  y  -  axis. 
Bearing  this  in  mind,  the  depression  along  the  ridges  at  tp  =  rr/  2.  3rr/  2  in 
figure  6a  and  subsequent  plots  is  understood  as  a  combination  of  Figures  5d  with 
the  shifted  versions  of  5b  -  5e.  The  configuration  of  scatterers  which  gave  rise  to 
Figure  6d  has  a  symmetry  line  along  <f>  -  rr/  4.  and  we  clearly  see  their  sym¬ 
metry  in  the  scattering  amplitude  plot.  The  sequence  of  Figures  6a  -  6d  is  a  pro¬ 
gression  toward  this  symmetric  arrangement. 

C.  Scattering  From  a  Plane. 

We  now  consider  adding  parallel  lines  of  scatterers  so  as  to  form  a  plane 
square  made  out  of  twenty  five  scatterers.  Figures  7a  -  7d  show  the  relevant 
sequence  of  plots.  It  helps  to  realize  that  the  final  square  configuration  (produc¬ 
ing  Figure  7d)  is  symmetric  around  =  0.  rr/  4  and  rr/  2  .  The  sequence  of  Fig¬ 
ures  7a  -  7d  shows  progression  toward  the  symmetric  result. 

D.  Scattering  From  a  Cavity. 

Finally,  we  consider  scattering  from  cavities  obtained  from  the  square 
configuration  of  Figure  7d  by  removing  scatterers  from  the  interior.  Comparing 
Figures  7d  and  8c  it  is  apparent  that  there  is  significantly  more  structure  in  the 
cross  section  for  scattering  from  a  cavity  than  for  scattering  from  the  rectan¬ 
gle.  The  increased  multiple  scattering  effects  in  the  cavity  enhance  the  subsidi¬ 
ary  maxima  at  0  =  rr/  4,  rr/  2  and  3rr/  2. 
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CALCUUTE  '$(0.0.0).  % (2.0.0),  '$(3.0.0)  BY  DIRECT  INTEGRATION 
ALL  OTHER  GREEN’S  FUNCTIONS  ARE  DETERMINED  BY  THESE  THREE 


P  a  l 


#‘(2p*l.  1.0)  ■  (.4-0)/ 2 
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f  a  ap-l 

m‘  *  3 


*8x(t+Zm+Z0)  »  -C/  (m'4-i) 
'S'im+ZC+Z,  0)  «  '8'(l-+Zm  +2,  o) 
m'  a  m'  +2 


'Sl(2p+ 2.0.0)  «  (*p.D+£-F)/ 

‘SM0.2JH-2.  0)  a  g'^^0.0) 

(gl(2P •►2.2.0)  *  D-'$l(2p+  2.0.0) 
'£*<2.2p-*-2.0)  »  <$x(2p+ZZ0) 
<$x(2p*Z  4.  0)  -  (0-  r$x(2p*2.  0.0) 

—4  <SM^»^2.0))/3 
«g»(4.  ^+2.  0)  a  ^‘(2p+2.  4.  0) 
r  -  2p 
m*  *  4 


#  ‘  (r  +2,  m-  +2.  0)  =  —C/  (m'+i) 

^5 1  (m'  +2.  f  +2.  0)  a  S'(M*'*2.0) 
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Figure  2. 


figure  5.  Polar  perspective  plot  of  I  f  (  0p  )  I  for  scattering 
scatterers  in  a  line  with  maximum  amplitude  l/m**! :  (a)  n 


Figure  5  (continued):  (c)  n  =  4,  l/maj  a  119  and  (d)  n  =  5,  i  I  ~  146. 
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e  7  (continued)  (c)  n  =  4  ,  I / mu  I  *  391;  (d)  n  =  5  .l/^l  =*  575  . 


Figure  0.  Polar  perspective  plots  of  i  f  (  ©,50  )  I  for  scattering  from  5x5  square  of 
point  scatterers.  In  (a)  the  center  scatterer  is  absent  and  I /ma, |l  ~  566.  In  (b) 
only  the  scatterers  on  the  perimeter  are  kept  and  I  f  I  «  445. 
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